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Abstract 

We propose a method to obtain the equilibrium distribution for positions and veloc- 
ities of a one-dimensional particle via time-averaging and Laplace transformations. 
We apply it to the case of a damped harmonic oscillator in contact with a thermal 
bath. The present method allows us to treat, among other cases, a Gaussian noise 
function exponentially correlated in time, e.g., Gaussian colored noise. We obtain 
the exact equilibrium solution and study some of its properties. 



PACS: 02.50.-r,05.40-a,05.40.Jc 



1 Introduction 



A Brownian particle is a massive particle (M), surrounded by a bath of lighter 
particles (m) [1]. Its motion is very irregular and varies rapidly with time. 
However, its averaged behavior can be well understood. The reason for this is 
that the mass ratio (m/M <C 1) leads to the separation of time-scales for the 
slow Brownian's degrees of freedom from the fast bath's degrees of freedom [2]. 
This allows us to construct adequate methods (e.g. Fokker-Planck equations) 
to describe the time evolution for the probability distribution function, or 
equivalently, for the momenta of the distribution [3,4]. The equilibrium dis- 
tribution is thus obtained as a function of the equilibrium momenta. We do 
not need to point out that the behavior of Brownian particles has been the 
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subject of many scientific studies since Einstein's seminal papers [5,6], already 
a hundred years old. Interesting extensions of this old problem have attracted 
the interest of many groups lately [7]. The treatment of non-standard forms 
of the noise function has also been successfully obtained recently [8] . 

In this paper, we calculate the equilibrium distribution for a damped har- 
monically bound Brownian particle as a direct integral, by means of Laplace 
transforms. However, a word is necessary to explain the use of the heavy ma- 
chinery of Laplace transforms upon such a well known model. 

In fact, the present model is only used as an illustration for the time-averaging 
method we propose. Applications for more interesting systems will follow. It 
basically consists in arranging the full dynamical behavior, expressed by the 
exact time-averaged distribution, into sums of integrals that can be proved 
either to vanish identically or to be easily calculable. Thus, the exact equi- 
librium distribution is obtained. Further extensions of the method are under 
way for non-Gaussian noise [9]. 

Technically, we define the equilibrium distribution in the spirit of Boltzmann's 
prescription [10] as a time-averaged distribution. The equilibrium distribution 
calculation is done by carrying out analytically the calculation for the time 
average of the distribution defined below, via the use of the stochastic functions 
x(t) and v(t), solutions of the stochastic dynamics. It reads 

1 T 

P eq (x, v) = lim - / dt <5(x- x(t))S(v - v{t)) >, (1) 





where <> means average over the noise rj(t). For ergodic systems fulfilling 
/(* -> oo) = lims/(s), 

s— >0 



where / is the Laplace transform of /, we can use an equivalent form for Eq. 1, 
namely 

oo 

P eq (x, v) = lim z [ dt e' zt < 5(x - x(t))S(v - v(t)) > . (2) 



As mentioned above, the functions x(t) and v(t) are solutions of the Langevin 
equations that follow from coupling the Brownian particle to a bath of lighter 
particles. We assume these to influence the Brownian dynamics via an effective 
non- white Gaussian noise force r](t). The effect of time-correlations leads to 
the renormalization of the oscillator frequency at equilibrium (all results are 
exact), discussed in the next section. 
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Our approach is straightforward and allows us to separate the reversible from 
the irreversible parts of the dynamics. The equilibrium distribution (the irre- 
versible part) is obtained by opening the delta-function on Eq. 2 by Laplace 
transformation into an infinite series, and by keeping only non- vanishing terms 
in it (in the limit z — > 0). As will be made clear in the following, the vanishing 
terms (in the limit t — > oo <^ z — > 0) do so due to their singular structure. 
They have poles (the same as those for the dispersion relation of the damped 
harmonic oscillator described by the Langevin Equation minus the random 
force) expressing the dissipation due to the damping. The surviving terms are 
the ones that do not carry the effect of the damping to the main propagator 
(as shall be defined below, proportional to z). The equilibrium distribution is 
obtained by summing over the non-vanishing terms. 

The method we propose could be summarized as follows: 

• Laplace transform of the Langevin equation and obtain x(z) and v(z); 

• Laplace transform of the Equilibrium distribution defined by Equation 1; 

• To eliminate the vanishing terms (in the limit z — > t — > oo); 

• To integrate over the the non-vanishing terms and obtain the equilibrium 
distribution. 

This method for obtaining the equilibrium distribution via time-averages, can 
be extended to different systems 3 . 

As stressed above, an important aspect of the method is that it allows us to 
study time-dependent processes (e.g., equilibration), since the effect of time 
comes from a single propagator. The long-time behavior of any dynamical 
functions for the system are easily singled out in the limit z — > <^ t — > oo. 

This paper is organized as follows. In Section 2 we define the model. In Sec- 
tions 3 we define, and Laplace transform, the time-averaged equilibrium dis- 
tribution. In Section 4 we discuss the irreversibility from the point of view 
of surviving and vanishing integrals. In Section 5 we select the non-vanishing 
contributions and obtain the equilibrium distribution in an analytical form. 
In section 6 we make some concluding remarks. 



3 Work is under way in order to extend the present results to more general forms 
of noise, in special for the case of non-Gaussian white noise. 
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2 Model 



2. 1 Langevin Equations and Laplace transforms 



We model a Brownian particle subject to a damped harmonic potential. The 
damping is proportional to the velocity. A Langevin force rj(t), Gaussian dis- 
tributed, effectively acts on the Brownian particle due to the interaction with 
the bath particles, not included in the dynamics. We assume that n(t) can be 
successfully Laplace transformed into fj(z). 

The system obeys the following coupled Langevin-like Equations 4 : 

mv(t) = -7«(t) - kx(t) + rj(t), (3) 
x(t) = v(t), (4) 

where r)(t) is a colored noise (time-correlated) distributed stochastic force. 
The initial conditions are given by 
x(t = 0) = x , v(t = 0) = v . 



Integrating Eqs. 3 and 4, we obtain that x(t) and v(t) can be expressed re- 
spectively as: 

t 

x(t) = x + J v{t')dt\ (5) 



and 



t t t' t 

mv (t) = mv - 7 / v{t')dt' - kx t - k [ [ v{t")dt"dt' + f rj(t')dt'. (6) 



o o 



The Laplace transform of fit) is defined as 

oo 

/(*) = Jdte- Zt f(t), (7) 
o 

4 Eq. 3 is strictly of Langevin type when the noise n(t) is white. A discussion on 
this can be found in reference [4], chapter IX. 
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with Re(z) > 0. 



By taking the Laplace transform of Equations 5 and 6, we obtain: 

x(z)=A(z)+B(z)fj(z), (8) 

and 

v(z) = C(z) + D(z)rj(z), (9) 



where 



a(z) = {z ; p)x « +v i do) 

z(z + B) + wf 



m ] 



B ( z ) = ZTTT^V-^i (11) 



z(z + (3) + wl 

ZVp - WqXq 

z{z + B)+w 2 

m~ l z 
z(z + B)+w 2 ' 



c(z) = j"°rf° xo . 2 , (12) 



D W = Z7TT^-^> (13) 



where wl = — and B — J -. 

Let's observe that only the coefficients A and C carry the contributions from 
the initial conditions. They are the memory terms: the vanishing of their 
contribution will make the process irreversible. 

The common denominator of the coefficients A, B, C and D is the damped 
harmonic oscillator's dispersion relation. Its poles are 



«± = -§ ±i\^4^7 2 , (14) 

where we assume that 4u;q — f3 2 > 0. This will be justified in the sequence, by 
the end of Section 5. 

We observe that the coefficients A and C carry information about the initial 
conditions, but not B and D, which are related to the stochastic part of the 
dynamics. The system will retain memory of its initial condition as far as 
terms depending on the coefficients A and C contribute to the final results. 
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2.2 Colored noise rj(t) 



We assume that the stochastic process r\ (t) is of the form of Gaussian colored 
noise [4], in which rj(t) is uncorrelated for t < 0. It is more general than white 
noise and reduces to it in the limit r — > 0. 

Some considerations have to be made about our choice of colored noise. For a 
closed system, noise has only internal causes. Using projection operator tech- 
niques [11] one obtains a generalized Langevin Equation [12] showing that the 
fluctuation-dissipation theorem holds near equilibrium and has to be consis- 
tent with a retarded dissipation of the form: 



v{t) = - J dt' 4>{t - t')v{t') + T](t), 



where 

(v(tHt + r)) 



0(t) 



knT 



However, we are assuming that the much simpler form of dissipation holds, 
namely Eq. 3. This can be justified by assuming that in a general situation (far 
from equilibrium), the fluctuation-dissipation theorem does not apply and the 
dissipation's memory kernel and the fluctuation properties of the noise need 
not to be related [13]. A method for solving the fully retarded dissipation for 
arbitrary noise can be found in reference [13]. 

Another interesting situation is the case of purely external noise. In this case, 
noise is not described by a closed Hamiltonian and dissipation and noise fluc- 
tuations are not related. 5 

Thus, we start from Eqs. 3 and 4 and study the long-time behavior of a 
non-Markovian coupled system (x, v). As for the noise properties, due to its 
Gaussian nature, only the second order cumulant survives (zero average dis- 
tributed): 

(V(t)) = 0; (15) 
(V(t)v(t)) = £e H *-''l/T. (16) 



5 We could imagine a simple system consisting of a charged particle immersed in a 
viscous fluid, driven by a strong external fluctuating electrical field. 
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As usual, for even / we have: 



i all pairwise 

(II viu)) = E (viuMu,)) ■ ■ ■ MvX*0>> 

i=1 arrangements 

for all possible combinations for I = 2p even (it will have combinations of 
pairs of r] functions). 

By Laplace transforming, the 77-average second cumulant reads: 



(fj(iqi + e)fj(iqj + e)) 



D 



[i(qi + qj) + 2e] [1 - riiq, + e)] [1 - r{iq 3 + e)] 



D 



T 



3 + r [i( qi + q 3 ) + 2e] - r 2 (ig, + e)(igr,- + e) 



2 [ [1 - r(ift + e)] [1 + r(i qi + e)] [1 - r{iq 3 + e)] [1 + r(i qj + e) 



.(17) 



Eq. 17 reduces to the white noise case 
D 

+ 9j) + 2e ' 



as t — > 0. 



The presence of time-correlated noise will affect the equilibrium distribution. 
As shall be seen, it has the effect of renormalizing the oscillator frequency. 
The white-noise equilibrium form is recovered when r — > 0. 



3 Equilibrium Distribution 



A convenient way of defining the long-time limit (t — > 00) for a convergent 
equilibrium distribution is to use 



P _ lim fS°dte-«f{t) 



for z > 0. This is equivalent to 
1 T 
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whenever this expression converges. In this case we write 

/(*-00) = /e, + A(i), 



with lim^oo A(t) = 6 . For z — ^, the average /f then reads 



eg 



lim 



ssyte-«f(t) 
6 /( 



00 c -zt 



dt 



lim 



f °°dte-t/ T f(t) 



T 



= f eq + lim 



lim 

<5^0 



/eg /< 



eg 



/eg 

< lim <5 



T^QO 

0. 



T 



The two forms of time averaging are equivalent. In the case the distribution 
defined on Equation 1 does not converge to an analytical form (as t — > 00), 
both definitions above will fail to give convergent values. 

Thus, we write the time-averaged equilibrium distribution as (see Appendix 
A): 



P eq (x,v) = \imz f dte- zt (5(x-x(t))5(v - v(t))} 



+00 



Umlim E , *QdP e ,^(-iQy(-jpy 

2 ^o^o,^„ J 2vr 27T U ml 



,m=0_* 



X 



n — ft dph 



f=l 271 h=l 27T z 



E/=i m + Eh=i Wh + + m)e 



1 



x (II £ (^/ + e ) \[v{Wh + t))- 
f=i h=i 



(19) 



The integration path for the q- and p- variables is also given in Figure B.l. Both 
upper and lower semi-hemispherical integration paths will contribute nothing 
to the total integral. We choose the upper path for our calculations. 



4 Origin of irreversibility at the limit z — > 



The equilibrium distribution we obtain corresponds to a sum of (infinite) prod- 
ucts of polynomial fractions with poles given by the dispersion relation, Equa- 

6 We assume that: V<5 > 0, 3M : Vt > M => \ A(t)\ < 5. 
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tion 14. One of those propagators, in Equation 19, stands out for being the 
only one that contains the effect of time: the main propagator G(z), given by 



G(z) 



z 



(20) 



z - [Kli + ■■■ + qi+Pi + ■■■ + Pm) + (I + m)e] ' 



In order to evaluate the equilibrium distribution (the analysis can also be 
extended to the non-equilibrium case z ^ 0), we need to study the poles and 
the limiting properties of the main propagator given by Equation 20. 

When integrating a given q or p-variable, by choosing the upper path on the 
complex plane shown in figure B.l, the pole associated with the propagator 
falls outside the integration path, contributing with no residues. However, after 
effecting all (m + I) variable calculations, values associated with the poles of 
the the rest of the integrand in Eq. 19 might be carried through into the 
denominator of G(z). It is essential for us to analyze G(z) at the limit when 



When the denominator in the main propagator is non-zero for z — > 0, it is 
clear that lim^o G(z) = 0. A non-zero result can only be obtained, in the limit 
z — > 0, when the sum of q and p-variable is completely eliminated from G(z). 
When all variables are eliminated, we have instead lim 2 ^ G(z) = lim^o f = 
1, and the integral yields a non-zero contribution. 

Since the factors A and C are not coupled with the stochastic force fj, the 
q and j9-variables associated with them will lead to residues only around the 
poles of the dispersion relation (— ik±). The pole's values will in the end be 
transported all the way to the denominator of the main propagator. So, by 
integrating over the poles of the functions A and C, a non-zero final value 
to the summation i [91 + ... + <#+ pi + ... + p m ] + (m + l)e, will give a zero 
contribution in the limit z — > 0. 

This is the simple way by which irreversibility manifests itself on the path to 
equilibrium: the memory of the initial conditions is lost as the inverse Laplace 
transform of terms such as 



z 



0. 



z 



z 



mpt 
2 



z — i(m (— m±j) 



z-m(3/2±i... 



e 



where m is an integer. 
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5 Important integrals 



The effect of poles from the dispersion relation in functions B and D can be 
eliminated since they are coupled with ^-functions, which are associated with 
a factor 8{iq + iq + 2e) dependency, coming from the stochastic properties 
of the ^-functions. When we integrate over the q and q variables, the 5(iq + 
iq + 2e) dependency eliminates the effect of these variables from the main 
propagator in the equilibrium limit z — > 0. Thus, only three possible couplings 
may contribute to the equilibrium distribution: BB, DD, and BD. 

In the following we list the non-vanishing contributions to the equilibrium 
distribution P eq (x,v). 



5. 1 Contribution from BB 



The typical contributing term for the x distribution is given by the BB con- 
tributing terms below (for details see Appendix B): 



-t-o 

/ 



t^ z-i( qi Z +q 3+ o) B{iqi + e)B ^ + + ^ + e)) = 

D X rS l^ + r 3 } {2l) 



z — io2[m + r 2 k — T7][m + r 2 k + T7] Icy k 7 

where o represents all the non integrated variables. 
Consequently, the only contributing terms will be even powers m = 29: 
m\ 1 Dm 2 T 2 D^y r 2 Dm t 3 D m 

'1-17™ 1 _2L\2 l OL ' nZ ' O Ji 2 ~~ 



2 f(^)! L [(m + r 2 A;) 2 -rV] 2A; 7 2k 2 1 

(20)! 1 Drrf T 2 D-f r 2 Dm t 3 D 6 

~ ~2M^ [(m + r 2 A;) 2 - r 2 7 2 ] ^1^7 ~ 2A; + 2 7 + ~2~^ 



5.2 Contribution from DD 



Similarly with the BB, the DD contributing terms are: 



/dp ■ dp ■ z 
^ z- i (p l+ p J+ o) D ^ + emP > + + + 6)) 
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'° r ""' r ' 2k +r] (22) 



2; — io 2[m + T 2 k — T7][m + T 2 k + T7] 7 7 



In the limit z — > 0, the contribution from the velocity momenta averages for 
the probability distribution is given by, for n = 2a: 

n\ 1 Dm r 2 Dk tD § _ 



2f(f)! L [(m + r 2 A;) 2 -rV] 2 7 2 7 2 

(2a)! 1 Dm r 2 Dk rD 

~ 2^P \{m + r 2 A;) 2 - r 2 7 2 ] ^ + ^7 2~ 



5.5 Contribution from BD 



Similarly with the above, we can calculate the cross contribution BD except 
that it has an odd power of the integrated variables (0(p)) in the numerator 
while the previous ones had even powers of it (0(1) and 0{p 2 ) respectively). 
Thus: 

/ ^ .-i^lp^o) 11 ^ + ^ + + + <)>, (23) 



the only important term 

+ F d Pj iz m~ 2 D Pj 

J 2n z — i(0 + o)(vj — in+.)(p-i — in-)(vi + iK+Mvi + 



-00 



is identically null because of the odd power of pj. So the equilibrium distribu- 
tion will be separated into a position term and a velocity term: 

P eq (x,v) = P eq (x)P eq (v). 
5.4 Equilibrium distribution 

Collecting all the no n- vanishing terms in the z — > limit, we write Eq. 19 as: 
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j-jQ) 20 m { 



i 



Dm 2 r 2 D^f T 2 Dm r 3 D 1 e 
29\ 2 e 9\ L [(m + r 2 k) 2 - r 2 7 2 ] [ ~2k^ ~ 2k + 2 7 + ~2~" * 



J-iP) 2 * (2a)! r 



1 



Dm r 2 L>A; rO 
2a! 2 a a! L [(m + r 2 k) 2 - r 2 7 2 ] ^ + 2~ 



TdQ Q 2 

J 2n eXVWX 2[(m + r 2 k) 2 -T 2 1 2 Y 2A; 7 2A; ' 2 7 



Dm 2 r 2 D-f T 2 Dm r 3 L> 

+ ^ + — ]>x 



r dP P 2 Dm r 2 Dk tD 

x y _ exp| ^ _ __ _ ^ _ ^ L __ + __ _j| 



Integrating the equation above, we finally obtain the properly normalized equi- 
librium distribution: 



[(m + r 2 /c) 2 — r 2 7 2 



\tt[ 



Cm 2 



T 2 D~f _|_ T ' 2 Dm 
fc 7 



+ r 3 D] 



m 2 [(l + r 2 ^ 2 ) 2 -r 2 /3V , 



&7 



+■ + T~ 3 -D] 



7 



X 



+ r 2 /c) 2 — r 2 7 2 ] 



7T|— + £ 

L 7 7 



- rD] 



( m 2 [(l + r 2 ^ 2 ) 2 -r 2 /3> , 

Dm r 2 £)l, _ nl eX Pi jDm | r 2 Dfc j J 



(25) 



On the limit r — > we obtain the distribution of displacement and velocity of 
a Brownian particle under a Gaussian white noise as expected: 



, A; 7 _22_kx^_ [rfvy 2 7 mv 2 
P eg (x,v) = \l—e-» 2 yj^e-n 2 



(26) 



Eq. 26 is the analytical time-averaged equilibrium distribution for the damped 
harmonically bound Brownian particle. 

The result from Eq. 25 has to be compared with exact results in the literature. 
The exact solution for a Langevin equation with time-correlated noise has 
been obtained by means of functional methods [14,15] or by path-integral 
approaches [16]. In order to compare our results with these references, we allow 
the system to be unconfined (we assume Eq. 25 holds for k — > and m — 1) 
and we integrate over the coordinate variable x. The reduced distribution 
obtained from Eq. 25 then reads 



7 [l + r 7 ] . 7 [l + T/3]t> 2 
exp{ } 
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This corresponds exactly (up to linear order on r) to the results in refer- 
ences [15,16] when we make the change D — > 2D due to slightly different 
definitions of the noise time-correlation. As for the result in reference [14], 
they are identical to all orders of r. 

A difference between our approach and these other authors' [14,15,16] comes 
from allowing another set of coupled variables (position {x}) to interact with 
our variable (velocity {v}). This leads to the appearance of terms propor- 
tional to kr 2 on the equilibrium distribution. As somewhat advanced in the 
Appendix A of reference [16], the equilibrium distribution only depends on 
simple combinations of small powers of r (on the exponent's argument). 



The fluctuation-dissipation ratio ^ plays the role of a thermal-bath temper- 
ature in the limit 7, D — > with = constant: 

7 rp _D 

2 7 



This is the reason we choose 4u; 2 > (3 2 in Section 2.1. 

However, for non-zero r the expression for the temperature grows a little more 
complicated. It will be given by 



, r r col-no l se_ £ t 1 ~ T P + ^j] 

B bath "2 7 [(l + T 2 u; 2 ) 2 -T 2 /3 2 r [ 8J 



From the form for the temperature and the equilibrium distribution, Eqs. 26 
and 28, we observe that the equilibrium frequency is renormalized to ui, given 
exactly by 



u 2_ u 2 [1 ~ T0 + r 2 ul} 
1_ [l-r 2 /5 2 + r^ 2 + r3^]- ^ 



As can be seen, in first order in r this corresponds to 



which reflects a decrease in the oscillating frequency due to damping. The time- 
correlations in rj(t) tend to make the velocities more persistent when they are 
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highest (near the equilibrium position) thus making the system getting further 
away from the origin and taking longer to come back. 

We are presently generalizing this method in order to treat other forms of 
potential gradients affecting the Brownian particle. 



6 Conclusions 

We present an approach for obtaining the equilibrium distribution of a simple 
system (in this case a Brownian particle under time correlated colored noise 
(with Gaussian white noise as a special case) via time-averaging (at the limit 
t — > oo) since taking the time-average is the most correct way of describing 
the long time equilibrium behavior of any physical system. 

Our model separates clearly the the irreversible contributions from the terms 
that keep track of the initial conditions (memory terms). It allows us to study 
each term contributing to the probability distribution for the position and 
velocity of the Brownian particle and to follow it as t — > oo. It is shown, 
by Laplace transforming the distributions, that the non-contributing terms 
(memory terms) disappear due to the vanishing of the main propagator as- 
sociated with them in the limit z — > <^ t — > oo. The non- vanishing terms 
can then be grouped together and integrated, yielding the exact equilibrium 
probability distribution, as has been confirmed by comparing with previous 
exact results in the literature [14,15,16]. 

In summary, the main advantage of the method is that it allows us to ob- 
tain the exact derivation of the equilibrium (and also the transient behavior 
if necessary, but not done in this paper) in a systematic way. It "opens" the 
problem into treatable parts (however ugly they look) that reach easily calcu- 
lable values at infinitely long times. 

Work is under way to extend the present results to other forms of noise. 
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A 



The detailed Laplace transformation of Equation 19 is done below: 

oo 

P eq (x,v) = \imz f e- zt (5(x-x(t))5(v - v(t)))dt 

z—*0 J 



= lim z I 

z^O J 



oo +00 
zt 



dte~ 



% 11 "0 ml 
lim, [die- I «±e* f ^ ± t*l ± ^ 



-00 
00 +00 
-zt 







/ =1 h=1 

00 +00 

limz / dte~ z ' 

z->0 



dQ iQ 


+ TdP 

x I e ipv 

J 2tt 

—00 


27T 


dQ &iQ 


+ r°dP 

x e iPv 




J 2tt 

—00 




- tia)S(t - t r 


dQ iQ 


+ r°dP 

x 1 ar JPv 


27T 


J 2tt 

—00 



E 



■(At)v m (t)) 



11 "0 m - 



f=i h=i 
HQY g (-tpy 



% 11 "0 mi 



x 



X 



+ ?° 1 rln m rlr, °r 1 °r m 

/nf n^/n*//n*- 

4,, f=i A7T h=i A7C J f=i J h =i 

e E a=1 (f-fi«)(%^)+Er=i( t - t ^)(^+ f ){ n x(t if ) n v(t mh )) 

f=i h=i 



v , dQdP eiQx+iPv {-iQl (-jP) 
z ^°i^=o_i 2vr 2tt /! m! 



/=i ~" fe=i 



2tt 



X 



2 — 



E l f =1 iqf + Y:T=iiPh + (l + m)e] /=1 



n%/+£) n^+^xA- 1 ) 



h=l 



The typical contribution for the equilibrium distribution for x is (integrations 
follow the path described in Figure B.l): 



dqi dqj 



I """" 

-oo 
+c 

/ 



— oc 
+oo 



2n 2n z — i(qi + qj + o) 
dqi dqj z 



B(iqi + e)B(iqj + e){fj(iqi + e)fj(iqj + e)) 



2tt 2tt z - i(q,i + ^ + o) (q { + + iK-){qj + + + <7j) + 2e] 
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m 



'-D 



z — z(0 + o) J 2tt [qj — z/c+][<7j — iK-](qj + in + )(qj + in J) 
D z 



2~jkz-i(Q + o) 



(B.l) 
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;. B.l. Integration path for the q or p-variables. 
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